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ABSTRACT 

We report the first Far Ultraviolet Spectroscopic Explorer [FUSE] measurements 
of diffuse O VI (AA 1032,1038) emission from the general diffuse interstellar medium 
outside of supernova remnants or superbubbles. We observed a 30" X 30" region 
of the sky centered at / = 315?0 and b — —41?3. From the observed inten- 
sities (2930±290(random)±410(systematic) and 1790±260(random)±250(systematic) 
photons cm^^ sr“^ in the 1032 and 1038 A emission lines, respectively), derived 
equations, and assumptions about the source location, we calculate the intrinsic inten- 
sity, electron density, thermal pressure, and emitting depth. The intensities are too 
large for the emission to originate solely in the Local Bubble. Thus, we conclude that 
the Galactic thick disk and lower halo also contribute. High velocity clouds are ruled 
out because there are none near the pointing direction. The calculated emitting depth 
is small, indicating that the O Vl-bearing gas fills a small volume. The observations can 
also be used to estimate the cooling rate of the hot interstellar medium and constrain 
models. The data also yield the first intensity measurement of the C II 3s ^Si /2 2p 
^P 3/2 emission line at 1037 A and place upper limits on the intensities of ultraviolet line 
emission from C I, C III, Si II, S III, S IV, S VI, and Fe TIL 
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1. Introduction 

In one popular description of the Milky Way’s interstellar medium, the Sun lies within a cluster 
of small, cool clouds (Lallement & Bertin 1992; Lallement et al. 1995; Linsky et af. 2000; Redfield 
h Linsky 2000), which are themselves embedded wdthin a large (approximately 70 pc in radius) 
bubble of highly ionized, rarefied gas (Cox h Reynolds 1987; Sno\vden et al. 1990; Warwick et al. 
1993; Frisch 1995; Snow^den et al. 1998; Sfeir, et al. 1999). The bubble is called the Local Bubble 
or the Local Hot Bubble and contains highly ionized gas, as evidenced by its | keV X-ray surface 
brightness and O VI absorption column density (Sanders et al. 1977; Fried et al, 1980; McCammon 
et al. 1983; Shelton & Cox 1994; Oegerle et al. 2001). Gas temperatures of 10^'^ to 10^ K have 
been calculated from collisional ionization equilibrium models of the observations. However, it is 
also possible that the gets has a lower temperature if it was heated by an ancient explosion and has 
subsequently adiabatically expanded and radiatively cooled (Breitschwerdt et al. 1996). 

Beyond the Local Bubble lies the Lockman layer of neutral hydrogen (Lockman et al. 1986), 
the Reynolds layer of ionized hydrogen (Reynolds 1993), and the Galaxy’s extended gaseous halo. 
The Local Bubble may be a bounded entity, though some recent evidence suggests a connection 
with the halo (Welsh et al. 1999). Like the Local Bubble, the lower parts of the halo or the upper 
parts of the thick disk are rich in highly-ionized gas. O VI, N V, and C IV ions have been found 
to extend to several kiloparsecs from the disk (Savage, Sembach, & Lu 1997; Savage et al. 2000; 
also Danly & Kuntz 1993) and X-ray shadowing observations show^ that some fraction of the ^ 
keV X-ray background originates beyond the Lockman layer (Snowden et al. 1991; Burrow^s & 
Mendenhall 1991, Snowden et al. 2000). We do not know if these ions reside in a hot plasma in 
collisional ionization equilibrium or if they are the cooled tracers of previous heating events. Such 
events may include halo supernovae, shocks due to fast infalling clouds, and galactic fountains. 

Much of w^hat is known about highly ionized gas in the Galactic interstellar medium (ISM) has 
been derived from observations of soft X-ray emission and absorption-line measurements of the O VI, 
N V, and C IV ions. Measurements of 0 VI emission can complement these data. Furthermore, the 
electron density and thermal pressure can be calculated from the ratio of the doublet intensity to 
the absorption column density, w^hile the plasma’s ionization history can be inferred by comparing 
the O VI intensity and high resolution soft X-ray spectra. Many investigators have searched for 
0 VI emission from the diffuse interstellar medium, but the only previous detection of Galactic O VI 
emission was made by the Hopkins Ultraviolet Telescope [HUT). In addition to observing superno\’u 
remnants and the Loop I superbubble, HUT observed the diffuse ISM toward I = 218'’, h = 
detecting the 0 VI doublet in the wing of the terrestrial Ly/? emission with a surface brightness 
of 19,000±5000 photons cm'"^ s"^ sr~^ (Dixon, Davidsen, & Ferguson 1996). Upper limits on the 
O VI surface brightness have been established for other directions (see Holberg 1986; Edelstein & 
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Bowyer 1993; Dixon, et al. 1996; Korpela, Bowyer, h Edelstein 1998), with most finding ~ 7000 
photons cm“^ sr~'^. Most recently, Edelstein et al, (1999) used Espectrografo Ultravioleta 
extremo para la observacion de la Radiacion Dijusa [EURD) data to find an upper limit of only 
1200 photons cm“"^ sr”^ with 90% confidence for a swath along the ecliptic, and Murthy et al. 
(2001) used low-resolution Voyager to find an upper limit of 2600 photons cm“^ s”^ sr“^ with 
90% confidence toward I = lll^, b = 51°. 

Here we report on the first observations of O VI (AA1032 and 1038) and C IT (AA 1037) 
emission from the diffuse ISM outside of supernova remnants or superbubbles made with the Far 
Ultraviolet Spectroscopic Explorer {FUSE). Dixon et al. (2001) have also used FUSE to observe 
diffuse 0 VI emission from a sight line in the northern Galactic hemisphere. FUSE has sufficient 
spectra] resolution to resolve easily the 0 VI emission lines from H I Lyman ^ and O I airglow lines 
several A away (see Sahnow et al. (2000) for a characterization of FUSF7s performance). These 
observations provide important constraints on models for the Local Hot Bubble and Galactic thick 
disk and halo. 

In section 2, we discuss the properties of the observed sight line, including the soft X-ray 
emission and dust opacity in this direction. In section 3, we describe the observations and 

data reduction techniques. Our measurements of the intensity, velocity, and width of the observed 
O VI and C II lines are presented in section 4. In section 5, we calculate the optical depth, intrinsic 
intensities, column density, electron density, thermal pressure, and emitting depth. The results 
are compared with information gained from other observations. We conclude that some of the 
observed photons probably originate in a thin region associated with the periphery of the Local 
Bubble, while the majority of the observed photons probably originate in the thick disk and lower 
halo. The path length of the emitting gas in the thick disk and lower halo is probably between 
ten and several hundred parsecs, and thus is very small when compared with the O VI scale height 
found by Savage et al (2000). Section 6 contains a summary of our conclusions. 


2* Observed Sight Line 

Our pointing direction is / = 315.00°, b — “41.33° (RA ” 22:44:30.0, Dec = -72:42:00.0). 
The observed intensity may include contributions from the Local Bubble and material out through 
the Galactic halo. Other than the Local Bubble, there are no known supernova remnants or 
superbubbles in this part of the sky. The pointing direction is near the edge of the Magellanic 
Stream (Wakker &: van Woerden 1997), but no emission is observed at the Stream’s velocity. 

Soft X-ray data provide some information on the conditions of the hot gas along the observed 
line of sight, although the plasma producing the soft X-rays is probably somew^hat hotter or more 
highly ionized than the plasma producing the O VI photons. Within a ~ 5° radius of the FUSE 
pointing, the average ROSAT \ keV surface brightness is 850 X lO"® counts arcmin"^ (Kuntz 
2000). Figure 1 shows the Snowden et al. (1997) \ keV map of the sky; the surface brightness at the 
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position of the FUSE pointing is similar to the average observed for latitudes b < -30^. The map 
also shows that the pointing direction is not toward an unusual portion of the southern sky. The 
observed region lies between extended regions of above- and below- average soft X-ray intensities, 
where the X-ray intensity gradient appears to be associated with absorption by the Galaxy’s H I 
layer (see Figure 1). 

IRAS maps of 100 //m emission (Wheelock et al. 1994) show that an irregular patch of far- 
infrared (FIR) cirrus, about 30' in diameter, overlaps our 30” X 30” FUSE field of view^ Atomic 
hydrogen maps (21 cm) yield Nhi ~ 4 xlO^^ cm“^, although, due to the coarse binning (Dickey 
h Lockman 1990; accessed through Sky View at http://skyview.gsfc.nasa.gov/skyview.html), sig- 
nificant variations are possible. Based on published photometry and spectral classification we have 
also derived estimates of the visual extinction in this region. Using a search radius of 60', centered 
on our pointing direction, we find 8 stars with spectral classification in the Michigan Spectroscopic 
Survey (Houk & Cowley 1975) and B-V measurements from Tycho photometry (Hog et al. 2000). 
Their spectroscopic and parallax information is listed in Table 1, Using spectroscopic parallaxes 
we find weak evidence for an obscuring “cloud” at ~200 pc causing a color excess of as much as 
E(B - V) ~ 0.07 magnitude, with no evidence for further reddening out to ^ 750 pc. This red- 
dening is consistent with the column density derived from the 21 cm data (Diplas & Savage 1994), 
and we assign it to the FIR (100 fim) cirrus patch. 

Far-ultraviolet emission (e.g., from 0 VI 1032 and 1038 A) originating from beyond the cirrus 
will be subject to scattering and absorption by the cirrus. However, because both the emission 
source and the cirrus are probably extended, many of the 1032 and 1038 A photons scattered 
out of the beam may be replaced by 1032 and 1038 A photons scattered into the beam. Because 
we do not know the details of either the source geometry or the FUV scattering characteristics 
in this direction, we can only estimate the effective amount of extinction experienced by photons 
originating from beyond the cirrus. In one limiting case, there are losses due to both scattering and 
absorption, so the average interstellar extinction curve can be applied to estimate the extinction near 
1032 and 1038 A. For E(B-V) < 0.07 and Ry=3.1 (the standard value for the diffuse ISM), we find 
that the extinction curve and parameterization presented in Fitzpatrick (1999) yields A 1035 < 0.99 
mag. Transforming this into an optical depth and applying the basic radiative transfer equation, 
this yiSd^an estimated loss of intensity of < 60%. In the other limiting case, all photons scattered 
out of the beam are replaced by photons scattered into the beam. In this case, there are no net 
losses due to scattering but there are losses due to absorption. For the cases in which scattering is 
inconsequential and absorption alone is very small or the sight line traverses a low density section 
of the cirrus, we take as the extreme lower limit a 0% loss in intensity. By combining these cases, 
we take as reasonable limits that between 40% and 100% of the 1032 and 1038 A intensity that 
enters the cirrus from beyond will emerge from the cloud and arrive at the observer. 
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3. Observations and Data Reduction 

The observations discussed here were acquired in 1999 September using the low resolution 
{LWRS; 30" x 30") aperture of the instrument. Summaries of the observatory and 

on-orbit characteristics can be found in Sahnow et al. (2000). The data sets used in this work 
are archived under the program identification 120509. This paper only presents data taken with 
detector 1 because detector 2 was off during much of the observations. After eliminating intervals 
when the spacecraft ^vas in the South Atlantic Anomaly and editing out event bursts (Sahnow et 
al. 2000), the exposure durations wdth the lA and IB detector segments were 219 and 233 ksec, 
respectively. Of these, 101 and 99 ksec, respectively, were acquired while the satellite was in the 
night portion of its orbit. The differences between day and night data are discussed below. 

The data were processed using version 1.7 of the standard CALFUSE pipeline (Sahnow et 
al. 2000; Oegerle, Murphy, & Kriss 2000), with four modifications: (1) the background noise 
originating in the detector wcts reduced by excluding counts which produced pulseheights outside 
the expected range for cosmic photons (the CALFUSE pipeline cutoffs were 4 and 12 in the standard 
arbitrary units), (2) the standard CALFUSE detector-background subtraction algorithm was not 
used, allowing us to calculate the uncertainties in our measurements more easily, (3) we extracted 
the spectra from the spatially corrected two-dimensional images of the detectors’ signal in units 
of counts reported as an intermediate CALFUSE pipeline data product, and in constructing the 
extraction windows, we used the the criterion that they must be extended far enough in the cross 
dispersion direction to at least encompass the entirety of each terrestrial airglow diffuse emission 
signal, (4) the LiF IB point source effective area curve, which corrects for the thin shadow^ of the 
repeller grid seen in point source spectra over the ~ 1155 to ~ 1170 A range (also known as the 
“worm”), was replaced with an effective area curve that does not include this correction. The 
relative wavelength scale of the resulting spectrum is accurate to five or six detector pixels, or 
about 0.035 A in the LiF lA spectrum. The absolute wavelength scale is less well known, and we 
discuss the determination of a wavelength (velocity) zero-point below. The intensity calibration 
was achieved using the best instrumental affective area estimates as of June 2000 (Sahnow et al. 
2000). The uncertainties in the solid angle of the aperture and effective area calibration are the 
largest sources of systematic error. The LWRS’s solid angle is thought to be accurately known to 
about 10%, while the effective area calibration is thought to be accurate to ~ 10% for LiF lA, SiC 
lA, and SiC IB, and ~ 20% for the LiF IB outside of the region of the “worm.” 

Figure 2 presents the spectra extracted from the detector 1 data in counts versus wwelength 
(i.e., no intensity calibration has been performed). LiF and SiC refer to the lithium fluoride and 
silicon carbide mirror and grating coatings used on the respective channels. The strongest emission 
lines seen in Figure 2 come from the Earth’s atmosphere. Across most of the bandpass, we see 
only noise due to the detector background, light scattered within the spectrograph, and cosmic 
background radiation. The small elevations in the intensity between ~ 1045 and ~ 1055 A in the 
LiF lA spectrum and between 1025 and ~ 1035 A in the SiC lA spectrum are due to light 
scattered within the spectrograph during the day portion of the orbit. This scattered light makes 
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it untenable to use the day time SiC lA data for analysis of the O VI emission. The night-only 
portion of the data provides less than half of the total exposure time. For this reason and because 
the SiC lA channel has a smaller effective area at 1032 A than does the LiF lA channel, the SiC 
lA data are not used in the O VI analysis. However, the night SiC IB data are used to search for 
S VI 933 and 944 A and C III 977 A emission since the LiF data do not cover these wavelengths. As 
can be seen in Figure 2, the LiF lA data show a peak in the vicinity of 1032 A that is potentially 
an O VI line, and another near 1037-1038 A that potentially contains a C II line and the other O VI 
line. 

We considered the possibility that some of the peaks observed in the 1030 to 1040 A region 
of the LiF lA spectrum could be due to light scattered within the spectrograph. Such scattering 
would be expected to distribute some light in the cross dispersion direction, as with the ~ 1045 to 
~ 1055 A contamination, which would be visible outside of the spectral extraction region. However, 
there is no excess of counts perpendicular to the dispersion direction in this region of the LiF lA 
spectrum. Thus, we rule out this type of scattering. 

We excluded the possibility that the observed 1032, 1037, and 1038 A emission lines could 
have come from the Earth’s atmosphere because of their continued strong emission in the satellite- 
night portion of the observations. Atmospheric airglow emission (such as the H I Lyman-/? and O I 
emission lines at rest w^avelengths of 1025.7, 1027.4, 1028.2, 1039.2, 1040.0, 1040.9 A) are roughly 
an order of magnitude stronger during the satellite-day portions of the observation than during the 
satellite-night portions. In contrast, the O VI and C II emission features in our data have constant 
intensity levels between the satellite-night and satellite-day portions of the observation. Thus, the 
emission features at 1032, 1037, and 1038 A are not due to terrestrial airglow’ (see Figure 3). 

We examined the possibility that scattered solar emission lines could have produced the fea- 
tures observed in the LiF lA data. The Sun’s ultraviolet spectrum has been recorded by the Solar 
Ultraviolet Measurements of Emitted Radiation [SUMER] instrument aboard the Solar and Helio- 
spheric Observatory [SOHO], The two strongest emission lines common to the SUMER bandpass 
and the combined FUSEIAY lA and LiF IB bandpass are the O VI 1031.9 A line and the 1175.7 A 
member of the C III triplet (Curdt et al. 1997). In the solar spectrum, the C III 1175.7 A line is 
0.86 times as bright as the O VI 1032 A line. Thus, if the 1032 A photons observed by FUSE czme 
from the Sun, then FUSE should have also observed 1175.7 A emission with 86% the intensity 
level of the 1032 A emission. But, the FUSE spectrum contains no emission features within at 
least 1 A of 1175.74 A. For A = 1175.74 ± 0.2 A, the intensity recorded by FUSE^ and measured 
with the techniques described below, is 360 ± 400(random) photons s”^ cm*“^ sr“^. This provides 
upper limits to the possible solar contamination of the LiP lA spectrum at 1032 and 1038 A of 
420 ± 470(random) and 210 d= 230(random) photons s“^ cm”^ sr”\ respectively. 

The LiF mirrors are physically situated on the shadow^ed side of the satellite, wiiile the SiC 
mirrors are on the sun-illuminated side. This motivated us to search for possible solar contamination 
in the SiC data. The C III 977 A emission line observed in the solar spectrum was also found in the 



-7- 


day time SiC IB spectrum but was not found in the night time SiC IB spectrum. The observed 
intensity is more than can be explained by terrestrial airglow. Thus, we conclude that the SiC 1 
data taken during the day portion of the satellite’s orbit ma}^ contain solar contamination, but the 
data taken at night contains no significant contamination. 


4. Spectral Analysis 
4.1. O VI Emission 

Table 2 presents the measured intensities and uncertainties, cr, the central velocities, the de- 
duced intrinsic full width at half maximum values, and the implied kinetic temperatures for the two 
O VI features. The values given in Table 2 were determined as follows. The background continuum 
was found by fitting a polynomial to the spectrum outside of the H I Ly/5 airglow, O I airglow, and 
next four strongest emission features (the 1032 A O VI emission line, the 1038 A O VI emission line, 
the 1037 A C II emission line, and a spurious ^1031 A feature). The continuum was subtracted 
from the spectrum. The intensities of the O VI emission features were measured by summing the 
signal over each line. The statistical (random) la uncertainty in the number of photons in each 
feature was calculated as the square root of the sum of the signal and background. The la sys- 
tematic uncertainty was calculated by adding in quadrature the 10% uncertainties in the aperture 
solid angle and flux calibration. The wavelength scale was corrected using an interpolation pinned 
to the rest wavelengths of the four strongest nearby airglow H I and O I emission lines (Morton 
1991, 2001) and the CALFUSE wavelengths at which these features were observed. The correc- 
tions applied to the 1032 and 1038 A emission features were —0.24 and -0.23 A, respectively. We 
then calculated the velocities relative to the geocenter and the Local Standard of Rest (LSR). The 
wavelength scale is probably accurate to within several detector pixels, or about 10 km s“^, and the 
centroid determination is probably accurate to within about 10 km as well. The line centroids 
were found by fitting Gaussian functions to the profiles. The observed widths are greater than the 
widths of the underlying cosmic spectrum because the FUSE instrument convolves the intrinsic 
spectrum with the instrumental point spread function and the large slit function. In an analysis of 
unblended O 1 and Ar I airglow lines in this and other data sets, we found that the instrumental 
effects are almost completely responsible for the the shapes and widths of the unblended airglow 
features (e.g., see the 0 I feature shown in Figure 4). Thus, we take the nearby 1039 A O I air- 
glow feature’s profile as the net instrumental function. We convolved this function with Gaussian 
functions of various widths, then compared the results with the observed 0 VI profiles in order to 
determine the intrinsic widths of the O VI emission lines. By equating the line dispersion with a 
thermal width we were able to calculate an upper limit on the kinetic temperature of the emitting 
plasma. 

The measured intensities of the O VI 1032 and 1038 A emission features are 2930±290 (random) ±4 10 (system 
and 1790±260(random) ±250 (systematic) photons cm~^ s“^ sr“^, corresponding to detections at 
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lOa and 7a levels of significance with respect to the random errors. Henceforth, the a associated 
with the random errors will be denoted by ^^(r)” and the a associated with the systematic errors 
will be denoted with by “(s)”. The velocity profiles of these lines are shown in Figure 4. For 
comparison, the (unresolved) 0 I airglow profile is also plotted. The 1032 and 1038 A features 
have similar central velocities, observed and intrinsic widths, and the observed velocity profiles 
have roughly similar shapes. The agreement corroborates our conclusion that the 1032 and 1038 A 
features are produced by a single species. 

The observed ratio of the intensities in the 1032 A and 1038 A lines (calculated from unrounded 
values; the values in Table 2 were rounded to 3 significant digits) is 1.64±0.29(r). The systematic 
uncertainties affect the 1032 and 1038 A emission line measurements in the same manner and thus 
the ratio of the 1032 and 1038 A intensities is not affected by the systematic uncertainties. In the 
limit of the optically thin case, as the optical depth approaches zero, the ratio should approach 
that of the statistical weights, 2. In the optically thick case, self absorption may change the ratio. 
The resulting ratio depends on the optical depth and the geometry of the O Vl-rich material. The 
ratio that we have observed lies 1.2 a below the theoretical optically thin value. The upper limits 
on the kinetic temperature surpass the temperature at which O VI is most abundant in collisional 
ionization equilibrium plasma, but the breadths of the O VI profiles could in large part be due to 
non-thermal processes such as turbulence or dispersion in bulk motions. 


4.2, C II Emission 

We have also observed emission from the 3s ^Si /2 2p ^F 3/2 transition of C II (A — 1037.02 A). 
The emission appears to be centered at a corrected LSR wavelength of 1037.16 A, The emission 
probably arises from the ^5x/2 to ^Ps/2 transition rather than the ^Sif 2 to ^P \/2 transition to 
ground state C II (A = 1036.34 A). This is probably because (1) the apparent wavelength is much 
closer to 1037.02 A than to 1036.34 A, (2) the A value for the ^5i/2 to ^Ps /2 transition is more 
than twice that of the ^5x^2 to ^Fx /2 transition, and (3) the ISM’s plentiful population of ground 
state C II should readily absorb 1036.34 A photons. 

Table 3 presents the observed intensity, cr, central LSR velocity, and width of the Gaussian fit 
to the 1037 A feature. The emission feature is redshifted relative to rest, but to a slightly lesser 
degree than the O VI lines. The C II feature is also significantly narrower than the O VI features, 
as might be expected from gas in a lower ionization stage. The observed FVVTIM is similar to the 
instrumental FWHM for diffuse emission that fills the large aperture, i.e., the C II emission line is 
poorly resolved. 
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4*3. Search for H 2 Fluorescence between 1031 and 1038 A: 

Hydrogen molecules in high rotational states have several possible transitions that would pro- 
duce 1032 to 1038 A photons. We have applied two tests to confirm that the emission lines observed 
by FUSE are indeed those of the O VI doublet and C II (^Si/2 ^^ 3 / 2 ) rather than fluorescing 

H 2 - In the first test, we set limits on the strength of the H 2 lines in the vicinity of 1032 A and 
1038 A by finding other lines in the bandpass that would be excited by the same mechanism 

and using the FF5F spectrum to measure their strengths. Morton &: Dinerstein (1976), Sternberg 
(1989), and Abgrall et al. (1993a, b) list the following H 2 emission lines: Lyman (6,0) P(3) tran- 
sition at 1031.19 A, Werner (1,1) P(2) transition at 1031.36 A, Werner (1,1) Q(3) transition at 
1031.87 A, Lyman (6,0) R(4) transition at 1032.35 A, Lyman (5,0) R(l) transition at 1037.15 A, 
and Lyman (5,0) P(l) transition at 1038.16 A. In order for an H 2 molecule to make the Lyman 
(6,0) P(3) transition, the ambient radiation field must have first pumped the molecule to the J=2 
rotational state, v' = 6 vibrational state, and Ylu electronic excitation state. The strength of 
the Lyman (6,3) P(3) (1172.38 A) transition arising from the same upper state should be 58% as 
strong as the Lyman (6,0) P(3) transition at 1031.19 A(Abgrall et al. 1993a). We have searched 
the wavelength-corrected FUSE spectrum for emission at this wavelength having the width of the 
nearest airglow line, finding an integrated intensity of -1300 ± 620 photons cm“^ s~^ sr”^. Thus, 
fluorescing H 2 molecules probably are not contributing 1031.19 A photons. We have applied this 
method to each of the H 2 transitions above. In every case, the signal is less than 2cr. We used 
the Sternberg (1989) model for our second test. In Sternberg’s spectrum of H 2 lines emitted by 
cool clouds bathed in ultraviolet light, the Q(l) Werner (0,0) line at ~ 1009.77 A and the P(3) 
Werner (0,1) line at ~ 1058.82 A should be at least 138% and 152% as bright as any of the H 2 
emission lines between 1030 and 1040 A. We have measured the intensities of these comparison 
lines, finding integrated intensities of 39 ±110 and 72 ±170 photons cm~^ sr“^, respectively. 
Thus, H 2 fluorescent emission is unlikely to be responsible for the intensities observed near 1032, 
1037, and 1038 A. 


4.4. Upper Limits on Other Emission Lines: 

We searched the LiF lA and IB spectra for interstellar emission from C I at 1122.260 A, Si II 
at 1023.700 A, S III at 1015.502 A, S IV at 1062.664 A, and Fe III at 1122.524 A and we searched 
the night time SiC IB spectra for interstellar emission from S VI at 933.378 and 944.523 A and 
C III at 977.020 A. No evidence for emission from any of these lines was found. The resulting 3cj 
upper limits are listed in Table 4. 
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5. Physical Characteristics 

In this section, we derive equations that relate the observed intensities of the O VI emission 
lines, the optical depths of the emitting regions in these lines, the intrinsic intensities of their 
emission, and the column densities of the lines. Additional equations relating the intrinsic emission 
intensity and the column densities to the electron density, thermal pressure and physical extent of 
the emitting region or regions are also presented. We then calculate the physical characteristics 
for the following cases of interest. In cases and #16, the above quantities are calculated from 
the observed line intensities and their ratio. In case #la, we assume that there is no extinction by 
the possible cirrus material. In case #16, we assume that the cirrus reduces the 1032 and 1038 A 
intensities by 60%. The uncertainties in the values calculated for cases #la and #16 are very large, 
primarily owing to the large uncertainties in the line ratio. Thus vve also consider cases #2, #3a, 
and #36, which draw the 0 VI column density from assumptions about the physical scenario rather 
than the line ratio. In case #2, the sole source of the observed emission is assumed to be the Local 
Bubble and therefore an estimate appropriate for the Local Bubble is taken for the column density. 
No cirrus extinction is assumed for case #2 because the cirrus material is thought to lie beyond 
the Local Bubble. In cases #3a and #36, the observed emission is assumed to originate partly in 
the Local Bubble and partly in the Galactic thick disk and lower halo. Therefore an appropriate 
estimate is taken for the column density. In case #3a, we assume that there is no extinction by 
the possible cTrrus material. In case #36, we assume that the cirrus reduces the 1032 and 1038 A 
intensities by 60%. Cases #2 and #3 are probably the most plausible. 

After presenting the results for each of these cases, we discuss five significant assumptions made 
in this analysis (especially the assumptions of negligible scattering of 1032 and 1038 A photons into 
the beam), consider their validity, and discuss the sensitivity of our results to these assumptions. 
Finally, we use the various derived physical characteristics to constrain the possible sources for the 
observed emission. 


5.1* Derivations 

Our first goal is to determine how the observed intensities in the 1032 and 1038 A lines are 
related to the intrinsic intensities. By intrinsic intensity, we mean the intensity of photons produced 
in the plasma. This intensity would have been observed if resonant scattering and extinction by 
dust had not changed the number of photons heading toward our detector. 

By using a self absorption analysis, we are able to link the observed and intrinsic intensities 
to one another via the optical depth term, and find a ready explanation for the difference between 
the observed line ratio and the optically thin value. Note, however, that the difference between the 
observed line ratio and the optically thin value is only marginally statisticalfy significant. In the 
self absorption scenario, photons emitted by the O VI ions may be scattered out of the beam by 
interactions with other 0 VI ions within the emitting region and residing along the line of sight, 
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while only negligible numbers of photons produced by other 0 VI ions in the emitting region are 
scattered into the beam. In section 5.3) w^e question this assumption and discuss how the calculated 
quantities would be affected if non-negligible numbers of 1032 and 1038 A photons were scattered 
into the beam. 


The standard radiative transfer equations governing emission produced within a slab of mate- 
rial for the case of negligible scattering into the beam and negligible extinction by dust are Spitzer s 
(1978) equations (3-1), (3-2), and (3-3) for the case where no externally produced radiation is in- 
cident upon the slab and is constant with respect to These equations yield the observed 

intensity per frequency interval for a given emission line: 


= r 

JO 




hi. 


(!-<!-"■) 


( 1 ) 


and the intrinsic intensity per frequency interval: 


hu = f ^drl - ( 2 ) 

Jo 

where u is the frequency, is the specific intensity per unit frequency, jV is the emissivity coefficient, 
ki, is the absorption coefficient, r' is the optical depth variable, and is the total optical depth 
through the emitting region. Following Spitzer, in our radiative transfer derivations is in units 
of ergs cm”^ 8"“^ sr“^ Hz"^, ji, is in units of ergs cm”^ s“^ sr”^ Hz”\ and ki, is in units of cm 
However, when presenting the results, we give the specific intensity in units of photons cm s 
sr-^ Hz-^ 

We combine Spitzer’s equations (3-16), (3-23), and (3-25) to yield an expression for The 
gas is assumed to be far from equivalent thermodynamic equilibrium though it may be in collisional 
ionizational equilibrium. Thus Spitzer’s ^ can be taken as zero, because the populations of the 
upper levels are negligible. Thus, 


= < <^{Ai/) >, 

TTleC 


( 3 ) 


where Nj is the column density of ions in the j state, Au = u - Vo, Uo is the frequency at line 
center, < <l>{Av) > is the line profile function (the path length-averaged normalized distribution 
with respect to frequency of the cross section per O VI ion), e is the charge of the electron, fjk 
is the oscillator strength for the chosen transition, is the mass of the electron, and c is the 
speed of light. From Spitzer’s equation (3-25), ^ = 2.654 X lO'^, where the units should be 


cm 

s 


If the particle velocity distribution is Maxwellian, then < <^(Ai/) >= ^ b ^ , Here b is the 

velocity-spread parameter, Ty is the FWHM velocity, and b = O-fiF^;. If w^e define Tq as the optical 
depth at line center, then 






rriecb 


= r, 


( \J\u \2 
^ b f . 


( 4 ) 
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We now integrate equations (1) and (2) with respect to v in order to determine and /{. 
Assuming that is constant over the frequency interval where 1 — differs from 0, 




3 Vo 


^Vc JQ 


Jo ^Vo Jo 


A At/ \2 

^ )d(Au). 


( 5 ) 


Equation (4) of Jenkins (1996) provides a solution to the integral in our equation (5) for r<, ~ 15: 
= /o°°[l - )]dx = ^ For x = we find 


j _ juc (-1)” j 


A u\y/u A 


(6) 


If there is dust between the emitting plasma and the observer, then equation (6) must be modified 
to account for the losses due to extinction. Equation (6) becomes: 


lo 


A 


( 7 ) 


where a is the ratio of the intensity that emerges at the near side of the dust to the intensity that 
impinges upon the far side of the dust. As discussed in section 2, if the O Vl-rich gas is beyond the 
IRAS ciTTus patch, then a is likely between 0.4 and 1.0. 


Assuming that Jvj^v is approximately constant over the range of v where ^ 0, the integrated 
intrinsic intensity is: 


^Vo Jo ^Vo Jo 


TnC 


^ ^ d{Au) 


To 


A 


( 8 ) 


We will use the ratio of equation (7) to equation (8) to calculate the intrinsic intensities from 
the observed intensities. This requires knowing the optical depth for each emission line: t"^(io 32) 
and 7*0(1038)- If column density is known, then equation (4) can be used to determine each 
Otherwise, the optical depths can be estimated from the ratio of the observed intensities in the 


1032 and 1038 A lines in the following manner. The oscillator strength for the 1032 A transition 
is twice that for the 1038 A transition and the ratio of is given by Spitzer’s equation (3-31), 

/ — Jt/o(io38) jg resulting 


7t/o(1032) 


thus, r„(io32) is equal to 2 r„(io 3 s) and A(io32)«..(io32) ' A(,o3*)«.o(io38) " ^a,,o32)' 

ratio of equation (7) for the 1032 A line to equation (7) for the 1038 A line is: 


A(1032) _ ^Q038)F'(^'''o(1038)) 
A(1038) '^(1032)F’('''o(1038)) 


( 9 ) 


This can be solved for Tp^josg). 

Next we determine the average electron density from the intrinsic intensity and the column 
density. The following equation motivated our earlier derivations of the intrinsic intensity and 
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column density. Shull & Slavin (1994) present a relation for the average electron density of an 
emitting region in which the kinetic temperature and emissivity per O VI ion are constant; 

, \ (/{5ooo) _ 

" {av), Novi ~ 

rj.x/2 /^1.392'\ 7,5000 

where is the electron-impact excitation-rate coefficient for the 1032 and 1038 A lines combined, 
Qovi{T) the Maxwellian-averaged collision strength for de-excitation of the O VI doublet, T is the 
gas temperature, the gas temperature in 10^ K, A^i 4 the column density in units of 10^"* cm“^, 
and /asooo the intrinsic O VI doublet intensity in units of 5000 photons cm' '2 s-i sr-^ Shull 
Slavin present a parameterization of the Maxwellian-averaged collision strength: Cloyj{T) = 
6.43[0.5584+ (0.1363y - 0.1158y2 + 0.279)^1 (y)e^ + 0.1158y], where y^EojkT and Ex{y) is the 
first exponential integral, 

The thermal pressure can also be calculated. Assuming a cosmic abundance of helium and 
complete ionization of both hydrogen and helium, 

Pth/k=^imrieT. (11) 

We expect the relevant temperatures to be near 3.2 X 10^ K because this is the temperature at 
which the fraction of oxygen atoms in the O VI ionization state is greatest (Shapiro & Moore 1976). 
The O VI ions could exist within plasmas of a range of temperatures. For example, the fraction 
of oxygen atoms in the O VI ionization state is at least 10% of its maximum value in collisional 
ionizational equilibrium plasmas for temperatures between 2.2 and 6.4 x 10^ K. The electron density 
calculated from equation (10) varies little across this range. 

The depth of the emitting region (A/ or J dl) is related to the column density through Nqvi ~ 
J novjdl^ where noyj is the density of O VI ions. Thus, A/ is approximately equal to Novil'^oyi- 
The volume density of O VI ions can be found from noyi ~ '^e ^ ^ ^ ^ > where rze, rifj, 

and no are the densities of electrons, hydrogen atoms, and oxygen atoms respectively. For a cosmic 
abundance of helium and complete ionization of both hydrogen and helium, is equal to 1.2nj/. 
Using the Grevesse and Anders (1989) solar abundance of oxygen, uq is equal to 8.5 X 10“^n7/. 
According to the Shapiro and Moore (1976) collisional ionizational equilibrium curve for oxygen, 
the maximum value of noyijno is 0.26. This occurs at T = 3.2 X 10^ K. The resulting equation 
for A/ is, 

A( = 5400 X (12) 

TIq 

where T is assumed to be 3.2 x 10^ K, Noyi is in cm”^, is in cm“^, and Al is in cm. Although 
most of the O VI ions probably reside within gas having nearly the maximum value of novifno^ 
some of the ions may reside in hotter or cooler gas. 


(0.00921 cm"^) 


^ovi{T) 
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5.2. Calculated Values for the Physical Characteristics 

In this subsection, we use equations (4), (7), (8), (9), (10), (11), and (12) to relate the ob- 
served intensities, optical depths, intrinsic intensities, colunan densities, electron densities, thermal 
pressures, and approximate emitting depth. The results are presented in Table 5 and Figure 5. In 
case #la, we take a = 1.0 and use the observed intensities and equation (9) to determine t"^(io 38)? 
finding 6i(r)* There is no systematic uncertainty in calculated from the ratio of 

the observed intensities. Then t^(io 32) is calculated from r^(io38)- We use the optical depths, ob- 
served intensities, and the ratio of equation (8) to equation (7) to calculate the intrinsic intensities 
for the emission lines. The resulting intrinsic intensity of the doublet is 7, 100^2^gQQ^jf^^^ ± lOOO(s) 
photons s”^ cm“^ sr”^. We calculate the velocity parameter from the velocity FWHM in Table 2. 
We use the velocity parameter and to calculate the column density in the ground state, finding 
N = 5.1^3‘^l^j X 10^^. Note that due to the uncertainty in the observed line ratio, the uncertainty 
in N is very large. We use the doublet Intensity, column density, and equation (10) to calculate the 
electron density for T = 3.2 x 10^ K, yielding = O.OOGO^q equation (10) and 

Figure 5 show, the electron density varies little between T = 2 x 10^ and 2 x 10^ K. We use equation 
(11) to calculate the thermal pressure at this temperature, finding Pthjk = 3, 700]|^j’?QQ|^j ± 520(s), 
Equation (121 is used to calculate the emitting depth if the gas is at this temperature, finding 
A/ = i 0.64(s) x 10^^ cm, or ± 21(s) pc. 

This process is repeated for a = 0.4 in case #li. Table 5 presents the results for each c^Lse. 
In the table, Al is presented in units of parsecs. The calculated uncertainties are rather large for 
some of the tabulated characteristics. The main source of the propagated uncertainties in cases 
#16 and 7^16 is the relatively large uncertainty in the line ratio. As a result, it is instructive to 
examine cases in which the column density (and hence the ratio of observed intensities) is set by 
the assumption of a particular geometry. In case the emission is assumed to originate in the 
Local Bubble. For this case we take the absorption column density to be Novi ~ 1.6 x 10^^ cm”^ 
(Shelton & Cox 1994). Only negligible extinction is expected, so a is taken to be 1.0. From the 
column density, we determine the optical depths and the ratio of intensities observed in the 1032 A 
and 1038 A emission lines. The total observed intensity is assumed to be conserved. W^e then 
calculate the other quantities using the same method as in cases #la and #16. The results are 
presented in Table 5 and Figure 5. 

In cases #3a and #36, the observed emission is assumed to originate in part in the Local 
Bubble and in part the Galactic thick disk and lower halo. For these cases, we take the absorbing 
column density to be Nqvj ~ 2.0 X 10^^ cm“^ from Savage et al. (2000). Because the previously 
discussed cirrus may lie either closer than the thick disk/halo emission region or beyond it, the 
physical characteristics are calculated using a = 1.0 (case #3a) and a = 0.4 (case #36). The 
calculations follow the same method as in case #2 and are presented in Table 5 and Figure 5. Note 
that in this model, the Local Bubble, thick disk, and lower halo are assumed to have the same 
thermal pressure. This is the pressure calculated from the equations and listed in Table 5. 
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5.3. Caveats 

In this analysis, we have made the following assumptions, which we will discuss in more detail 
in the remainder of this section. 1.) The contribution of 1032 and 1038 A photons emitted by other 
O VI ions in the “cloud” but outside of the beam and scattered into the bearri via interactions with 
O VI ions is assumed to be negligible when compared with the number of such photons scattered 
out of the beam via such interactions. 2.) The contribution due to scattering of background 
ultraviolet continuum photons is assumed to be negligible. 3.) Losses due to other types of particles 
interspersed with the O VI ions are assumed to be negligible (though the calculations do consider an 
external region of extincting particles). 4.) The emissivity per O VI ion is assumed to be constant. 
5.) Furthermore, by using equations (10) and (11), we make the assumption that the temperature 
and noyijno ratio are near their values at the peak of the ion fraction curve constructed by Shapiro 
and Moore using collision al ionization equilibrium. 

If the first assumption is not valid then the intrinsic doublet intensities will be less than 
the calculated values. The degree of scattering into the beam is not known. If the number of 
photons produced by the O VI “cloud” and scattered into the beam via interactions with 0 VI ions 
compensates exactly for the number produced by the O VI “cloud” and scattered out via interactions 
with O VI ions, then the intrinsic intensities would not follow equation (7), but would be equal 
to the observed intensities after correction for dust extinction. The results of such a scenario are 
within the error bars on our calculated doublet intensities in each of our cases (i.e. #la, #16, #2, 
#3a, and #36). If the number of photons produced by the O VI “cloud” and scattered into the 
beam due to interactions with O VI ions were to be even larger than the number produced by the 
O VI “cloud” and scattered out, then the true intrinsic intensity w^ould be less than the observed 
intensity after correction for dust extinction. In this scenario, we would expect the 1032 A intensity 
to be at least twice as strong as the 1038 A intensity because O VI ions scatter 1032 A photons more 
readily. Such a scenario is marginally excluded because the observed ratio of the 1032 and 1038 A 
intensities (1.64 ± 0.29(r)) is at least 1.24 sigma from the theoretical value (>2). In cases #la 
and #16, if some of the photons scattered out of the beam have been compensated for by photons 
scattered into the beam, then the true optical depths and column density would be greater than the 
calculated values. For cases #2, #3a, and #36, the column densities were assumed and the optical 
depths were found from them. If some of the light scattered out of the beam were replaced by 
light scattered into the beam, then the ratio of the observed to the intrinsic intensity would be less 
than that calculated from equations (7) and (8) using the calculated from the assumed column 
density. In all cases, if the intensity scattered into the beam is not negligible when compared with 
the intensity scattered out of the beam, then the true electron density and thermal pressure would 
be less than the calculated values. 

It is difficult to evaluate the second assumption, primarily because the geometric relationship 
between the stellar continuum sources and O VI ions is not well knowm. The following is an order 
of magnitude estimate. We assume that the O VI ions on the line of sight are bathed in the 
incident ultraviolet continuum field, whose strength is: = 4.38 x 10”^^ ergs cm^^ sr”^ 
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Hz"^ (Mathis, Mezger, & Panagia 1983). The volumetric intensity absorbed by the 0 VI photons 
will be Iai> = ’ The ions will quickly re-emit this volumetric intensity, 

but a fraction of the re-radiated photons will arrive at the observer w^hile other photons will be 
scattered out of the beam due to absorption by intervening O VI ions or wall be absorbed by other 
intervening particles. The observed intensity [Isco) can be found by solving Spitzer’s equation (3-1) 
or our equation (1), where is substituted for ji, and w^here an extinction factor, a, is included. 
The result is 

Isco = ^h<.F{To). (13) 

We use equation (13) to create rough estimates of the possible contribution of scattered ultra- 
violet continuum photons to the observed intensity. The results are presented in the last column 
of Table 5. The systematic uncertainties in Isco are bound to be large. As a result, we are not 
prepared to make confident corrections to the observed 0 VI intensities at this time. However, it is 
logical to recognize that if some scattered ultraviolet continuum photons are being observed, then 
the true intrinsic intensity, electron density, thermal pressure, and emitting depth are less than the 
values presented in Table 5. In particular, if the estimated scattering is used, then the doublet 
intensities, electron densities, and thermal pressures in the most plausible cases, cases ^^2 and #36, 
would decrease from the quoted values by 4.0% and 17%, respectively. The emission depths would 
increase by 4.0% and 20%, respectively. 

In the third assumption, there are no other absorbing particles interspersed with the O VI ions. 
If this assumption is false and if the calculated extinction by the cirrus is not an overestimate, then 
the effective extinction w^ould be greater than assumed. Thus, the intrinsic intensities, electron 
density, and thermal pressure would be greater than the calculated values. The emitting region 
would be smaller than the calculated value. 

In the fourth assumption, the emissivity per O VI ion is assumed to be constant. In case 
#3 where some of the emission is from the Local Bubble and some is from the halo, the thermal 
pressure varies along the path length. Thus the electron density and emissivity vary along the path 
length. The higher pressure gas in the Local Bubble would be more emissive per ion than the low^er 
pressure gas in the thick disk and halo. As a result, in case #3, the pressure listed in Table 5 would 
lie between the pressure of the gas in the Local Bubble and the pressure of the gas in the thick disk 
and halo. This logic also applies to case #1. 

If the assumptions regarding the temperature and ionization state of the gas are invalid, then 
the electron density, thermal pressure, and emitting depth are affected- If the gas is near the 
specified temperatures but out of collisional ionizational equilibrium, then the emitting depths 
would be larger than calculated from equation (12). Alternatively, if the fraction of oxygen ions 
in the O VI stage is near the peak in the Shapiro and Moore curve, but the gas temperature is 
much low^er (or somewhat higher) than assumed, then the emissivity would be lowrer (or higher) 
than assumed and so the calculated electron density would be higher (or lower) than assumed. 
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5.4. Possible Locations of the Emitting Gas 

In this subsection, we critique possible emission sources. We first examine the possibility that 
the observed emission originates entirely in the Local Bubble. ROSAT observations of soft X-ray 
emission from the Local Bubble suggest that the thermal pressure is ~ 15, 000 K cm”^ (Snowden 
et al. 1998). This value was found by fitting the soft X-ray spectrum with model spectra for solar 
abundance, collisional ionization equilibrium plasma and is appropriate for the ~ 10^ K, X-ray 
emitting plasma. Cool clouds within the Local Bubble have been shown to have lower thermal 
pressures (i.e. Jenkins (1998)). We assume that the difference between the thermal pressures of the 
hot and cool gas is balanced by other forms of pressure, such as magnetic or turbulent pressure. 
Hence, the thermal pressure in the cool gas does not constrain our analysis. The pressure of the 
soft X-ray emitting gas is within the range of pressures found for case #2. 

If we assume that the O Vl-bearing gas is in pressure balance with the X-ray emitting gas and 
that the 0 VI ions predominately exist within gas of temperatures between 2.2 and 6.4 x lO^K, then 
the electron density is about x 92x(2 2 toT4xio^K) “ 0-012 to 0.036 cm”^. Using equation (10) with 
these temperatures and electron densities and the Shelton & Cox (1994) value for the Local Bubble 
column yields an intrinsic doublet intensity between ~ 440 and ~ 1300 photons cm“^ s~^ sr“^. Not 
only is this range substantially lower than the intrinsic doublet intensity calculated for case #2, 
it is well below the observed doublet intensity! Scattered ultraviolet continuum photons probably 
cannot account for the intensity discrepancy because the relatively small Local Bubble O VI column 
density is ineffective at scattering the assumed incident photon field. Thus, we conclude that the 
Local Bubble, which is assumed to be in pressure balance and collisional ionizational equilibrium, 
cannot account for the bulk of the observed O VI photons. An additional source is needed. 

The photons which do originate in the Local Bubble probably come from a thin zone of hot 
gas. Using a Local Bubble pressure of 15,000 K cm“^ and equations (10), (11) and (12), we 
calculate the thickness of the O Vl-bearing gas as a function of temperature. If the temperature is 
r = 3.2 X 10^ K, then the fraction of oxygen atoms in the O VI ionizational state is at its maximum 
and Al = 1. 1 pc. If the temperature is T = 2.4 x 10^ K or T* = 6.4 x 10^ K, where the abundance 
of O VI ions is 10% of its maximum value, Al = 7.9 pc or 23 pc, respectively. Thus it is possible 
that the O VI ions reside in a thin shell on the edge of the Local Bubble. Such a structure is 
consistent with the lack of O VI column density towards stars located within the interior of the 
bubble (Oegerle et al. 2001). 

If observed emission originates in the Local Bubble, thick disk, and halo, then the derived 
pressure should be between the thermal pressures of the hot gases in these regions. As mentioned 
above, the thermal pressure of the hot gas in the Local Bubble is thought to be about ~ 15,000 
K cm“^. The thermal pressure in the thick disk and halo ranges from 0 at large distances from 
the Galactic midplane to ~ 13, 000 K cm”^ within 2 kpc of the plane in regions where the total 
pressure is predominately in the form of thermal pressure (Boulares & Cox 1990). The derived 
pressures for cases #3a and #3b are within the allowed range. 
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From Table 5, depending on dust extinction, the thickness of the emission zone would be 12 to 
29 pc if r = 3.2 X 10^ K. If the temperature is between 2.4 and 6.4 x 10^ K, then the values would 
be larger by a factor of up to 7 and 20, respectively. On a Galactic scale, these values are not large. 
It is possible that the O VI in the thick disk and lower halo resides in midscale structures scattered 
about the thick disk and lower halo. Such a scenario would be consistent with the large scale height 
(2.7 kpc) and high degree of inhomogeneity found in the Savage et al. (2000) mini-survey. 

We now examine the expected velocity for the pointing direction. In general, the velocities of 
the highly-ionized material in this part of the sky can be understood by assuming that a thickened 
disk of highly-ionized material roughly corotates with the underlying thin disk of low-ionization 
gas (e.g., Lu et al. 1994; Sembach &: Savage 1992). Our measured O VI velocities differ from the 
velocities expected from co-rotating gas wdthin 10 kpc of the plane by over 100 km sec~^. This 
discrepancy does not rule out the possibility that most of the emitting material resides within the 
thick disk / halo. It is quite plausible that the 0 VI resides within shock-heated gas traveling over 
100 km sec“^ relative to the speed expected from galactic rotation. 

A third hypothetical location for the observed material is high- or intermediate-velocity clouds 
and their environments, but this hypothesis appears to be contradicted by the 21 cm observations. 
Wakker et al. (2001) have presented Parkes H I 21 cm observations in the direction of the star 
HD 3175, which lies ~ 14^ from the FUSE pointing, showing an H I intermediate-velocity cloud 
ulsr +73 km s ^ with A^(H I) 10^^ cm However, the H I spectra from the southern 
H I 21 cm emission survey of the Institute Argentine de Radioastronomia (Arnal et al. 2000; see 
also Morras et al. 2000) show no evidence for an H I intermediate- velocity cloud over the velocity 
range ulsr +50 to +90 km s ^ within ~ 1?5 of of the FUSE pointing (R. Morras 2000, private 
communication). 


6. Summary 

We have observed emission by O VI ions in the Galaxy’s general diffuse ISM outside of super- 
nova remnants or superbubbles. We have ruled out the possibilities that the observed photons are 
stellar, instrumental or terrestrial airglow\ The observed intensities, velocities, and intrinsic widths 
are presented in Table 2. We have used the observed properties to estimate the column density, 
intrinsic intensity, average electron density, thermal pressure, and depth of the emitting gas. The 
results are presented in Table 5. 

The doublet intensity is larger than expected for emission by the (assumed quiescent) Local 
Bubble alone. Thus, we surmise that some of the observed intensity must arise in the Galactic 
thick disk or lower halo. There are no H I clouds in the southern H I 21 cm emission survey of 
the Institute Argentine de Radioastronomia near this direction at similar velocities. Thus, high 
velocity clouds are probably not the source of the observed emission. The emitting material appears 
to be confined to relatively shallow regions such as the edge of the Local Bubble or 10 to 600 pc 
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thick structures in the Galactic halo. 
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Fig. 1. — Top: The ROSAT | keV surface brightness map of the southern Galactic pole reproduced 
from Snowden et al. (1997) , courtesy of Steve Snowden. The Galactic latitude is -90® at the center 
of the map and increases outward. Lines of latitude are drawn every 15°. The Galactic longitude is 
0° at the top of the map. Longitude increases counterclockwise, with meridians drawn every 30°. 
The units on the sidebar are 10”^ counts s~^ arcmin”^. The dark arc running across the lower 
right quadrant indicates a region lacking sufficient ROSA T All Sky Survey data. The ~ keV surface 
brightness for I = 315°, I = -41° is approximately average for the region of the southern Galactic 
hemisphere aw^ay from the plane. Bottom Left: Enlargement of the upper right corner of the X-ray 
map. A circle is drawn around / = 315°, b = —41°. Bottom Right: DIRBE-corrected IRAS 100 
lim map, reproduced from Snowden et ah 1997. 
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Fig. 2. — The counts spectra from the I — 315^, h = “41.3^ pointing. Each spectrum is logorithmi- 
cally scaled in counts, includes satellite day and night data, and is binned by 11 pixels (~ 0.075 A 
for the LiF lA). The strong, narrow peaks are terrestrial airglow emission lines. The elevations in 
the LiF lA count rate at short wavelengths and the SiC lA countrate at long wavelengths are due 
to elevated scattered light levels. The small elevations in the intensity between ~1045 and ~1055 A 
in the LiF lA spectrum and between ~1025 and ~1035 A in the SiC lA spectrum are due to a 
localized band of scattered light cutting across the spectrum. 
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Fig. 3. — The LiF lA spectrum, binned over 11 pixels (corresponding to 0.075 A or about one fifth 
of the width of the large aperture) and plotted relative to the uncorrected CALFUSE wavelengths. 
Top: spectrum from the satellite night portion of the data. Middle: spectrum from the satellite 
day portion of the data. Bottom: spectrum from the day and night data. The emission lines 
associated with the Earth’s atmosphere (each denoted by 0) are quite strong in the total and 
day-only spectra, but much weaker in the night-only spectra. The ~1032 and ^1038 A emission 
features of interstellar O VI have similar brightnesses in the night-only and the day-only spectra. 
The 1037 A feature of interstellar C II is also present in all spectra. The 1031 A feature appearing in 
the total and day-only spectra is absent in the night-only spectra and is not statistically significant 
in the day-only data, thus is probably not cosmic. 
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Fig. 4. — The 1032 and 1038 A 0 VI and 1039 A O I profiles versus LSR velocity and binned by 
5 detector pixels. Top: The intensity in the O VI 1032 A emission feature. For plotting purposes, 
the curve hcts been shifted upwards by 65 photons cm“^ sr”^ (km Middle: Intensity in 

the O VI 1038 A feature, shifted upwards by 30 photons s“^ cm”^ sr“^ (km s“^)~L The 1032 A 
and 1038 A features have similar central velocities, widths, and roughly similar envelopes. Bottom: 
To illustrate that the instruinent optics convolve the intrinsic profile with a top hat function, the 
observed airglow O I 1039 A profile, scaled by 0,3 and shifted to the right by 65 km s“^, is also 
shown. 
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Fig. 5. — The electron density as a function of temperature for an O VI plasma, assuming values 
calculated directly from the data (case #la, shaded region), values calculated from the data and an 
assumed O VI Local Bubble column density (case ^2^ dashed line), and values calculated from the 
data and an assumed Local Bubble plus thick disk and lower halo column density (case #3a, thick 
solid line). Overplotted are lines of constant thermal pressure {Pth/k = 1.92neT, dotted lines). 
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Table 1: Stars within 60 arcmin of / = 315.0°, b = —41.3° 


Star 

Ang. Sep. 

Sp. Class 

Phot. Dist. 

Par. Dist. 

E(B-V) 


(arcmin) 


(pc) 

(pc) 


HD214691 

37.0 

G5V 

40 

49ti 

0.05 

HD2 13591 

51.0 

GOV 

53 

47^2^ 

-0.02 

HD213928 

53.0 

F7/F8V 

65 

77is^ 

0.01 

HD214778 

40.0 

G2V 

97 

249tei^ 

-0.05 

HD214522 

59.0 

F8/G0V 

134 


-0.02 

HD2 13742 

41.0 

KlIII 

167 

1751?? 

0.07 

HD216379 

58.0 

K2III 

381 

600 _92 

0.07 

HD215970 

39.0 

G8/K0III 

741 


0.02 



Table 2. 0 VI Observations 


Characteristic 

0 VI 1032 A feature 

0 VI 1038 A feature 

Intensity (photons s”^ cm*“^ sr“^) 

2930 

1790 

Intensity (photons s“^ cm”^ sr“^) 

2930 

1790 

a due to random uncertainties (photons s“^ cm“^ sr“^) 

290 

260 

a due to systematic uncertainties (photons s~^ cm”^ sr”^) 

410 

250 

Center LSR Velocity (km s”^) 

+70 

+60 

Intrinsic FWHM (km s“^) 

~ 110 

~ 110 

Kinetic Temperature (K) 

< 4 X 10® 

rsj 

< 4 X 10® 

r\j 
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Table 3: Characteristics of the C II Emission Feature 


Characteristic C II 1037 A feature 

Intensity (photons s“^ cm“^ sr”^) 1990 

a due to random uncertainties (photons s~^ cm“^ sr“^) 260 

a due to systematic uncertainties (photons s“^ cm”^ sr“^) 280 

Center LSR Velocity (km s“^) +40 
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Table 4: 3a Upper Limits for Other Transitions 


Species 

Rest Wavelength 

(A) 

3a Upper Limit^ 
(photons cm”^ sr“^) 

S VI 

933.378 

3600 

S VI 

944.523 

3600 

C III 

977.020 

2600 

S III 

1015.502 

780 

Si II 

1023.700 

390 

S IV 

1062.664 

350 

C I 

1122.260 

580 

Fe III 

1122.524 

580 


Note - a: As described in the text, these measurements are also subject to the systematic uncer- 

tianties of 14%. 
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